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Based  on  Fourier’s  law  and  the  Seebeck  effect,  this  paper  presents  a  mathematical  model  of  a  Thermoelec¬ 
tric  Generator  (TEG)  device  using  the  exhaust  gas  of  vehicles  as  heat  source.  The  model  simulates  the 
impact  of  relevant  factors,  including  vehicles  exhaust  mass  flow  rate,  temperature  and  mass  flow  rate 
of  different  types  of  cooling  fluid,  convection  heat  transfer  coefficient,  height  of  PN  couple,  the  ratio  of 
external  resistance  to  internal  resistance  of  the  circuit  on  the  output  power  and  efficiency.  The  results 
show  that  the  output  power  and  efficiency  increase  significantly  by  changing  the  convection  heat  transfer 
coefficient  of  the  high-temperature-side  than  that  of  low-temperature-side.  The  results  also  show  that 
with  variation  in  the  height  of  the  PN  couple,  the  output  power  occur  a  peak  value,  and  the  peak  value 
decreases  when  decreasing  the  thermal  conductivity  of  the  PN  couple,  and  increases  when  increasing 
the  Seebeck  coefficient  and  electric  conductivity  of  the  material.  Meanwhile,  a  maximum  output  power 
and  efficiency  of  a  TEG  appear  when  external  resistance  is  greater  than  internal  resistance.  This  is  differ¬ 
ent  from  a  common  circuit,  and  with  the  increment  of  ZT,  the  maximum  value  moves  toward  the  direc¬ 
tion  of  an  increasing  ratio  of  external  resistance  to  internal  resistance.  Finally,  we  propose  a  new  idea  to 
reform  our  experiment  design  to  achieve  better  performance. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Transportation  vehicles  powered  by  combustion  engines  are  one 
of  the  significant  causes  for  the  petro  fuel  crisis  and  intensifying 
environment  pollution.  From  the  perspective  of  the  energy  balance 
of  a  combustion  engine,  the  power  efficiency  accounts  for  only 
30-45%  of  the  total  amount  of  heat  generated  by  the  fuel  used. 
The  rest  of  the  energy,  except  for  a  10%  or  less  loss  to  overcome  fric¬ 
tion,  is  not  fully  used.  It  is  mainly  emitted  to  the  atmosphere 
through  heat  dissipation  and  circuit  cooling  exhaust.  During  this 
process,  40%  of  the  energy  is  lost  through  gas  emission  waste  [1], 
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Capturing  a  small  portion  of  the  available  energy  with  thermoelec¬ 
tric  devices  can  reduce  engine  loads  and/or  alternator  size  thus 
decreasing  pollutant  emissions  and  fuel  consumption  [2],  Thermo¬ 
electric  devices  can  also  be  used  in  conjunction  with  traditional 
C02  reduction  approaches  (reducing  friction,  decreasing  pumping 
losses,  and  improving  combustion  efficiency)  to  further  reduce 
greenhouse  gas  emissions  [3].  Previous  studies  had  estimated  that 
if  the  amount  of  driving  energy  required  for  the  production  of 
electrical  power  (approximately  6%),  could  be  recovered,  this  would 
result  in  an  overall  reduction  of  fuel  consumption  by  10%  [4],  It 
makes  sense  that  recycling  and  reusing  waste  exhaust  gas  can  not 
only  enhance  fuel  energy  use  efficiency,  but  also  reduce  air  pollu¬ 
tion.  Thermal  power  technology  such  as  the  TEG  arouses,  therefore, 
significant  attention  worldwide.  TEG  is  a  technology  for  directly 
converting  thermal  energy  into  electrical  energy.  It  has  no  moving 
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Nomenclature 

Abbreviation 

h 

convective  heat  transfer  coefficient  (W/(m2  K)) 

TEG 

Thermoelectric  Generator 

C 

specific  heat  capacity  (kJ/(kg°C)) 

MPPT 

maximum  power  point  tracking 

M 

mass  flow  rate  (g/s) 

P 

output  power  (W) 

Symbols 

row 

n 

efficiency  (%) 

m 

line 

Subscript 

n 

total  number  of  PN  couples 

f 

exhaust 

T 

temperature  (°C) 

c 

cooling  fluid 

Q 

quantity  of  heat  (W) 

1 

hot  side 

S 

Seebeck  coefficient  (pV/K) 

2 

cold  side 

/l 

thermal  conductivity  (W/(mI<)) 

P 

P  material 

A 

cross  sectional  area  (m2) 

n 

N  material 

s 

cross  sectional  area  of  one  kind  of  material  in  one  unit 

i 

in 

(m2) 

0 

out 

a 

electrical  conductivity  (s/m) 

y 

y  direction 

H 

height  of  either  P  or  N  thermoelectric  materials  (m) 

overall 

total 

P 

electric  resistivity  (Clm) 

K 

thermal  conductance  (W/K) 

Superscript 

Ri 

internal  resistance  (£2) 

J 

j'th  unit 

r2 

external  resistance  (£1) 

parts,  is  compact,  quiet,  highly  reliable  and  environmentally 
friendly.  Because  of  these  merits,  it  is  presently  becoming  a  notice¬ 
able  research  direction  [5],  Several  years  ago,  some  car  manufactur¬ 
ers,  such  as  HI-Z  in  the  USA,  Nissan  in  Japan,  BMW  in  Germany,  had 
successfully  developed  a  thermoelectric  engine  to  recycle  vehicle 
exhaust  gas.  However,  the  quantity  of  heat  recycled  is  quite  limited ; 
in  general  less  than  5%.  This  is  somewhat  limited  by  thermoelectric 
materials.  Although  challenges  exist,  it  may  be  possible  to  recover 
significant  amounts  of  exhaust  waste  heat  with  thermoelectric  de¬ 
vices,  particularly  as  research  of  materials  and  application  of  heat 
transfer  technologies  continue  to  improve  the  performance  of  these 
elements.  Continued  research  in  this  area  is  further  encouraged  as 
China’s  current  national  goal  and  the  US  Department  of  Energy’s 
commitment  to  reducing  America’s  dependence  on  foreign  oil 
gathers  more  attention  [6], 

With  the  rapid  development  of  semiconductor  technology  and 
the  growing  crisis  of  decreasing  fossil  fuel,  more  and  more  attention 
has  been  drawn  to  this  area.  But,  what  we  must  emphasize  is  that 
heat  recovery  efficiencies  are  still  low  now,  both  on  a  system  basis 
(<2.5%)  and  on  a  thermoelectric  generator  basis  (<3.5%),  despite 
great  progress  made  over  the  last  decade  in  the  thermoelectric 
materials  figure  of  ZT  merit  values  [7],  The  figure  of  merit  value  is 
a  non-dimensional  number  that  provides  an  overall  assessment  of 
the  thermo-element’s  electrical  conversion  efficiency.  Although 
material  selection  and  ZT  values  are  important  to  the  performance 
of  the  thermoelectric  generator,  power  output  also  greatly  depends 
on  the  availability  of  heat  in  order  to  maintain  a  large  temperature 
difference  between  the  hot  and  cold  side  of  the  materials  element 
[7],  Bell  [8]  pointed  out  two  important  pathways  that  would  lead 
to  additional  applications  of  thermoelectric  (TE)  devices.  One  is  to 
promote  the  intrinsic  efficiencies  of  TE  materials.  From  the  stand¬ 
point  of  current  new  thermoelectric  materials  research  and  devel¬ 
opment,  a  large  number  of  studies  concerning  thermoelectricity 
focus  on  how  to  improve  heat-to-electricity  conversion  efficiency 
of  thermoelectric  materials;  that  is,  how  to  improve  the  thermoelec¬ 
tric  ZT  figure  of  merit  of  these  materials.  Nanotechnology  [9,10], 
novel  technology  in  solid  state  physics  and  semiconductor  physics 
[1 1,12]  is  employed  to  explore  this  exciting  field.  Another  approach 
is  to  improve  the  way  existing  TEGs  are  currently  used.  In  addition 


to  the  research  of  thermoelectric  materials,  reasonable  thermal  de¬ 
sign  and  management  of  the  thermoelectric  generator  is  equally 
important  for  improving  generating  performance  [13],  The  key  fac¬ 
tor  should  be  the  usage  of  economical  and  efficient  of  heat  sources. 
From  the  perspective  of  system  analysis  and  optimization,  Rowe 
and  Gao  [14]  provided  a  practical  procedure  for  optimizing  TEG 
module  geometry  guided  by  the  “economic  factor”.  Gou  et  al.  [5] 
tested  results,  and  discussions  showed  the  promising  potential  of 
using  a  thermoelectric  generator  for  low-temperature  waste  heat 
recovery,  especially  in  the  industrial  field.  More  than  2%  efficiency 
is  achieved,  increasing  slightly  with  a  boost  of  temperature  differ¬ 
ence.  Hsu  et  al.  [15]  proposed  a  concept  of  “effective  Seebeck  coef¬ 
ficient”,  which  discussed  the  inconsistency  between  the  theoretical 
Seebeck  coefficient  and  the  measured  one.  Applying  an  appropriate 
pressure  on  a  TEG  module  was  suggested  to  improve  performance. 
Bismuth  telluride  material  was  used  in  their  experiment  to  verify 
their  conclusion.  Liang  et  al.  [16]  also  constructed  an  analytical 
model  of  a  parallel  TEG  by  theoretical  analysis,  to  investigate  how 
contact  effect  reduces  the  output  power  of  a  parallel  TEG  and,  fur¬ 
ther,  to  verify  the  model  by  experimental  testing.  Hsiao  et  al.  [17] 
had  constructed  a  one-dimensional  thermal  resistance  model,  and 
the  results  were  verified  by  experiments.  The  model  results  showed 
that  a  TEG  module  presented  better  performance  on  the  exhaust 
pipe  than  on  radiator.  Kim  [18]  had  also  proposed  an  experimental 
method  to  extract  the  intrinsic  and  extrinsic  Seebeck  coefficients 
and  resistances  of  a  TEG.  Crane  and  Bell  [19]  had  studied  the  influ¬ 
ence  by  changing  the  temperature  of  exhaust  gas  and  the  mass  flow 
rate  experimentally,  and  gave  an  optimizing  method.  Eakburanawat 
and  Boonyaroonate  [20]  described  a  thermoelectric  battery  charger 
powered  by  thermoelectric  (TE)  power  modules.  It  might  not  be 
economical  initially,  but  the  results  derived  here  may  lay  a  founda¬ 
tion  for  further  investigations  of  TE  battery  charger  systems.  Yu  and 
Chau  [1  ]  had  combined  MPPT  and  Cuk  converter  technology  with  a 
TEG  device,  and  achieved  20%  improvement  in  performance.  Espin¬ 
osa  [21  ]  had  proposed  an  idea  that  connected  different  PN  materials 
that  had  different  “best  working  temperature  ranges”  in  order  to 
realize  the  effect  of  temperature  gradient  optimization,  and  to 
get  a  better  performance.  He  also  tested  the  thermoelectric 
performance  of  many  kinds  of  materials’.  Chen  et  al.  [22,23]  also 
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explored  the  performance  of  the  TEG  with  multi-element  thermo¬ 
electric  equipment  and  built  a  model  of  a  two-stage  semiconductor 
thermoelectric-generator  with  external  heat-transfer.  Multi-stage 
models  consist  of  different  temperatures,  such  as  low,  medium, 
and  high.  Thermoelectric  modules  were  also  proposed  by  other 
scholars,  such  as  Xiao  et  al.  [13,24]  and  Chen  et  al.  [23],  Their  results 
showed  that  reasonable  thermal  design  of  a  thermoelectric  module 
could  take  full  advantage  of  the  characteristics  of  each  thermoelec¬ 
tric  material  and  effectively  improved  the  performance  of  power 
generation. 

Several  authors  had  demonstrated  the  use  of  thermoelectric 
elements  for  automobile  heat  recovery,  installing  the  devices  in 
the  vehicle  exhaust  [25-28],  Hsu  et  al.  [29,30],  Masahide  et  al. 
[31],  Yodovard  [32],  Yu  and  Chau  [33],  Xu  et  al.  [34],  Gou  et  al. 

[5],  and  Suter  et  al.  [35]  applied  thermoelectric  generation  technol¬ 
ogy  to  environmentally-friendly  vehicles.  Both  waste  heats  from 
vehicle  exhaust  and  engine  were  converted  into  electrical  power 
to  provide  auxiliary  power  for  the  vehicle.  To  achieve  better 
performance  or  efficiency,  many  scholars  designed  various  lands 
of  integrated  systems,  such  as  combining  with  photovoltaic 
[36,37]  or  Organic  Rankine  Cycles  (ORCs)  [38],  Many  scholars  also 
studied  thermoelectric  cooling  technology  to  cool  important  vehi¬ 
cle  devices,  or  for  comfort  [39-42],  Most  existing  researches  are 
experiments  with  low-and-medium  temperature  heat  sources,  or 
analyzing  the  influence  of  external  conditions  such  as  temperature 
and  mass  flow  rate  of  both  cold  and  heat  sources.  However,  theory 
research  and  internal  optimization  of  TEG  devices  and  other 
related  contents  is  far  less.  It  is  noticed  that  the  exhaust  gas 
temperature  is  not  constant,  in  fact,  reducing  along  with  the  TEG. 
The  temperature  profile  inside  the  TEG  module  and  its  changing 
behavior  is  important  for  designing  a  thermoelectric  device.  If 
the  inlet  temperature  profile  of  the  device  can  be  achieved,  it  can 
provide  great  help  to  optimize  the  device  and  to  analyze  the  influ¬ 
ence  of  various  factors.  To  this  end,  a  simple  thermodynamic  mod¬ 
el  is  proposed  that  can  be  used  to  calculate  the  temperature  profile 
in  the  TEG  module,  and  to  theoretically  predict  output  power  and 
efficiency.  The  effects  of  various  factors  on  the  performance  of  a 
TEG  also  have  been  analyzed. 

2.  Mathematical  model 

Thermoelectric  Generator  (TEG)  device  consists  of  many  small 
couples  of  P  type  and  N  type  thermoelectric  materials  assembled 
according  to  definite  permutations  and  combinations.  Completely 
symmetrical  device  fluctuation  is  assumed,  regardless  of  parts 
studied,  as  shown  in  Fig.  1.  We  start  our  study  from  a  single  PN 
unit  considered  to  have  the  same  length  scale.  The  entire  TEG, 
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therefore,  can  be  divided  into  many  small  units.  The  outlet  temper¬ 
atures  of  the  former  unit  are  treated  as  the  inlet  temperatures  for 
the  next  unit;  thus  the  calculation  can  be  done  continuously.  Some 
assumptions  are  made  as  follows: 

(1)  Steady  state  in  thermal  and  flow. 

(2)  The  thermal  resistance  of  copper  chip  is  ignored,  as  it  is  quite 
small. 

(3)  The  gas  in  the  gaps  is  ignored. 

(4)  No  contact  thermal  resistance. 

(5)  Heat  conduction  is  perpendicular  to  flow  direction,  the  heat 
conduction  in  the  axial  direction  is  omitted,  and  the  heat 
radiation  is  also  omitted. 

(6)  The  thermal  conductivity  of  the  PN  unit  is  constant. 

(7)  Thomason  heat  is  ignored. 

(8)  It  has  r  rows  and  each  row  has  m  couples,  so  the  total 
number  of  PN  couples,  n,  can  be  defined  as  n  =  r  ■  m.  PN  cou¬ 
ples  are  connected  in  series.  We  assume  the  materials  of 
each  row  are  exactly  the  same,  so  performance  of  each 
row  can  be  the  same.  We  can  simplify  the  three-dimensional 
model  into  a  two-dimension  model. 

As  shown  in  Fig.  1,  the  TEG  works  between  a  heat  source  (ex¬ 
hausted  gas)  and  cold  source  (ambient  air).  The  heat  is  trans¬ 
ferred  from  the  heat  source  at  temperature  of  Tf,  and  released 
at  temperature  of  Tc.  In  this  process,  part  of  the  heat  is  con¬ 
verted  into  electricity  that  is  transmitted  in  a  form  of  output 
power  or  electromotive  force. 

In  the  present  model,  the  Seebeck  effect  and  Fourier’s  law  of 
heat  conduction  are  assumed  to  be  dominated  in  a  representative 
size  of  unit  as  shown  in  Fig.  1.  In  thej  unit,  the  energy  conservation 
equations,  including  P  and  N  type  materials,  can  be  given. 
Assuming  T\  be  the  averaged  temperature  of  the  hot  side  of  one 
unit  length  in  the  axial  direction,  the  heat  transfer  rates  (in  y 
direction)  through  each  branch  of  PN  type  thermoelectric 
materials  andqj,,  at  the  root  side,  respectively,  can  be  expressed 

§l,  W 

^  =  -^A--'LA%\yo  (2) 

Let  T{  be  the  averaged  temperature  of  the  cold  side  of  one  unit 
at  the  outer  end  of  y  =  H\  the  heat  transfer  rates  through  the  PN 
thermoelectric  materials  qj,2  and  q’n2,  respectively,  can  be  expressed 


Fig.  1.  A  schematic  diagram  of  the  theoretic 
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(3) 

<4  =  -  4/n  h  (4) 

where  S'  is  the  Seebeck  coefficient,  /  is  the  output  current,  )!y  is  the 
thermal  conductivity  of  the  y  direction  of  P  or  N  type  thermoelectric 
materials;  Ap(An )  is  the  cross  sectional  area  of  a  unit  of  P  or  N  type 
thermoelectric  materials;  dT/dy  is  the  temperature  gradient;  H  is 
the  height  of  either  P  or  N  thermoelectric  materials.  The  relation¬ 
ship  between  the  temperature  gradient  with  joule  heat  per  unit 
height  in  PN  couples  is  given  by: 


■  d2r_/VP 

pAdy2~  A? 

-LA^M 

ny  dy2  K 


where  p/  is  the  electric  resistivity  of  the  material. 
The  boundary  conditions  are  given  by: 


(5) 


(6) 


y  =  0;  T=TJ1  (7) 

y  =  H-  T  =  Tj2  (8) 


From  (3-6),  we  derive: 


•  dr  i2Pp(y  -  h/2)  4A(r2-n) 
V*P  -  Ap  +  H 

dT  I2p>n(y  -  H/2)  4/n(n-n) 

'^Andy~  An  +  H 


(9) 


(10) 


Eq.  (9)  together  with  Eqs.  (1 ),  (2),  and  (10)  together  with  Eq.  (3),  (4), 
the  unit  absorbed  heat  q[(y  =  0)  and  unit  released  heat  q'2  (y  =  H)  in 
j  unit  can  be  obtained: 


<7i  =  (Hi  -  ^)IT\  +  K(T\  ~  Ti)  ~  /2r,i/2  (H) 

Q2  =  (4  -  4)^1 + k(j{  -  Ti ) + /2Ri/2  (i2) 


As  the  temperature  of  fluid  must  be  changed  after  heat  transfer, 
therefore,  the  following  equations  can  be  obtained: 


<li  =  Mi  | 

(17) 

fTi  ti0  +  e\  ccMc(t{0-tQ 

(18) 

where  C  is  the  specific  heat  capacity,  and  M  is  the  mass  flow  rate. 
tj  indicates  the  temperature  of  the  exhaust  gas  and  t{  is  the 
temperature  of  the  cooling  air. 

First,  we  assumed  all  units  are  independent  and  then  divided 
the  load  resistance  into  n  portions.  Each  portion  is  allotted  to  each 
PN  couple,  so  we  can  calculate: 

(19) 

<*+<>  <*+<) 

where  R2  is  load  resistance  (external  electric  resistance),  n  is  the  to¬ 
tal  number  of  PN  couples.  It  is  seen  that  electric  current  I  is  still  a 
function  of  T{,T{.  Eqs.  (15)— (19)  constitute  a  complete  problem  to 
be  solved  with  four  unknowns  of  t}-0.  T\ ,  V2.  Vm,  as  tj,,  tJd  is  the  initial 
value,  or  the  value  of  last  unit.  However,  an  iteration  calculation  is 
needed.  Accordingly,  we  can  calculate  the  electric  current  /,  and  all 
the  other  unknown  temperatures  for  each  unit.  As  we  know,  the 
electric  current  I  should  be  equal  in  a  serial  circuit;  however,  the 
obtained  electric  current  /  for  each  PN  couple  might  be  different. 
To  resolve  this  issue,  we  sum  up  the  currents  of  all  portions,  and 
then  take  an  average.  The  averaged  value  is  put  back  into  the  calcu¬ 
lation  until  a  tolerance  is  satisfied. 

The  output  power  of  each  PN  couple  can  be  expressed  as: 


F=(4-S'1)/T' 

-(Sp2-Sn2)lV2-l2R\ 

(20) 

and  its  efficiency  c; 

V=P'/q'1 

in  be  given  by: 

(21) 

where 


4a,  4a,  piH  pLh 

H  +  H  ’  1  Ap  +  A„ 

K’  is  the  thermal  conductance  of  j  unit;  R',  is  the  electric  resistance 
of  j  unit  of  PN  couple.Note  that  the  absorbed  heat  at  the  hot  side  q', 
from  exhaust  gas  can  also  be  written  as: 


q,1=hiA'1 


(13) 


The  overall  power  output  POVerait  and  efficiency  of  power 
generation  rjovemU  can  be  obtained: 


Poveral,  =  (22) 

»  (23) 


3.  Results  and  discussion 


Similarly,  the  released  heat  tf2  at  the  cold  side  to  ambient  air  can 
also  be  written  as: 

cf2  =  fi2A2  (V-i,  -  ^  (14) 

where  A, .  AJ2  are  the  heat  exchange  areas  of  the  hot  side  and  cold 
side  in  the  j  unit,  respectively,  hi  is  the  convective  heat  transfer 
coefficient  of  the  hot  side;  h2  is  convective  heat  transfer  coefficient 
of  the  cold  side. 

So  the  following  equations  can  be  obtained: 


Mj| 

(4^-n) 

=  {%,-^)‘A+k(t{-V2) 

-I2R{/2 

(15) 

Mil 

l  =  (4-4)/Ti+K(ri -Ti) 

-I2^/  2 

(16) 

Fig.  2a  shows  how  the  mass  flow  rate  of  exhaust  gas  varies  for 
the  federal  test  procedure  (FTP-75)  city  driving  cycle.  Mass  flow 
rate  ranges  from  5  g/s  up  to  45  g/s.  Fig.  2b  shows  the  variation  in 
exhaust  gas  temperatures  exiting  the  vehicle  catalytic  converter 
for  the  FTP-75  drive  cycle  [17],  The  temperatures  fluctuate  around 
500  °C,  but  do  not  fluctuate  as  much  as  the  mass  flow  rates. 

Table  1  lists  the  parameters  for  PN  materials  used  in  the  present 
paper.  The  other  parameters  such  as  the  temperature  of  the 
exhaust  gas  and  convection  heat  transfer  coefficient  are  listed  in 
Table  2.  In  other  subsequent  analysis  without  special  description, 
the  parameters  are  taken  from  Tables  1  and  2. 

Fig.  3  shows  the  calculation  result  of  temperature  profiles  inside 
the  TEG  device  with  the  parameters  listed  in  Tables  1  and  2.  The 
temperature  of  exhaust  gas  decreases  from  500  °C  to  383.7  °C; 
cooling-air  temperature  raises  from  30  °C  to  86.6  °C;  temperature 
of  semiconductor  material  ranges  from  294.6  to  254.2  °C  on  the 
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Table  1 

Parameters  of  PN  materials. 


Parameters 

P  type 

N  type 

Seebeck  coefficient  $ 

40  +  0.20T  (flV/K) 

40  +  0.20T  (pV/K) 

Electrical  resistivity  p1 

1.04  x  10~5  (fim) 

1.04  x  10~5  (Om) 

Thermal  conductivity  f! 

1.5  (W/(mK)) 

2.5  (W/(mK)) 

Height  H 

0.005  (m) 

0.005  (m) 

Sectional  area  AP(A„) 

0.01  x  0.01  (m2) 

0.01  x  0.01  (m2) 

Table  2 

Other  parameters. 

Exhaust  inlet  temperature  Tf 

500  (°C) 

Cooling-air  temperature  Tc 

30  (“O 

Convection  heat  transfer  coefficient  of  hot  side  h, 

100  (W/(m2  K)) 

Convection  heat  transfer  coefficient  of  cold  side  h2 

100  (W/(m2  K)) 

Exhaust  mass  flow  rate  Mf 

20  (g/s) 

Cooling  air  mass  flow  rate  of  i 

Dne  side  Mc 

20  (g/s) 

Ratio  of  external  resistance  to 
Each  unit’s  area  A 

internal  resistance  R2/R i 

0.015  x  0.025  (m2) 

Total  area  Aovera„ 

0.15  x  0.5  (m2) 

Note:  R,  =  . 


Number  of  PN  couple  (j) 


Fig.  3.  Variation  of  temperatures  of  each  unit  for  air-cooling. 

hot  side,  while  from  226.3  to  215.4  °C  on  the  cold  side;  output 
power  is  65.6  W;  efficiency  is  2.66%.  We  can  conclude  the 
utilization  rate  of  waste  heat  of  exhaust  gas  is  quite  low  and 
directly  manifested  in  the  small  temperature  difference  between 


inlet  and  outlet  exhaust  gas  of  the  TEG.  Although  the  efficiency 
of  the  thermoelectric  is  general  low  at  the  present  stage,  lower  out¬ 
put  power  and  efficiency  are  mainly  caused  by  a  low  temperature 
difference  between  the  hot  side  and  cold  side  of  the  thermal  mate¬ 
rial  (which  is  almost  1  \1  compared  to  the  temperature  difference 
between  the  exhaust  gas  and  ambient  air,  as  Fig.  3  shows.  So, 
improving  the  coefficient  of  convective  heat  transfer  of  fluid  from 
a  heat  transfer  perspective,  and  developing  new  methods  of 
thermoelectric  transition  from  a  thermodynamic  perspective  are 
quite  necessary. 

3.1.  Effects  of  temperatures  and  mass  flow  rates  of  cooling-air 

The  variations  of  performance  with  temperatures  and  mass  flow 
rates  are  shown  in  Fig.  4.  Both  the  output  power  and  efficiency 
decrease  approximately  linearly  with  the  increment  of  cooling 
temperature.  Because  the  rise  in  cooling  temperature  makes  the 
temperature  difference  between  the  two  sides  of  PN  couple 
decrease,  according  to  the  Seebeck  effect,  resulting  in  the  decre¬ 
ment  of  voltage,  leading  to  the  reduction  of  output  power,  we  thus 
conclude  that  at  lower  cooling  temperatures,  the  better  the  perfor¬ 
mance.  With  an  increment  of  cooling-air  mass  flow  rate,  both  the 
output  power  and  efficiency  increase  nonlinearly.  When  the  mass 
flow  rate  exceeds  40  g/s,  the  improvements  of  output  power  and 
efficiency  are  both  quite  small,  while  bigger  mass  flow  rates  always 
cost  more  power.  There  must  be  an  optimal  cooling  mass  flow  rate 
value  according  to  all  parameters. 

3.2.  Effect  of  mass  flow  rates  of  exhaust  gas 

The  variation  of  performance  with  exhaust  gas  mass  flow  rates 
is  shown  in  Fig.  5.  With  an  increment  of  exhaust  gas  mass  flow 
rate,  the  output  power  and  efficiency  only  increase  noticeably 
when  the  mass  flow  rate  is  low;  while  the  waste  heat  absorption 
of  the  TEG  is  insufficient  when  mass  flow  rate  is  high.  Further, 
when  the  exhaust  gas  mass  flow  rate  increases  from  lOg/s  to 
30  g/s,  output  power  increases  only  20  W.  If  we  use  three  of  the 
same  devices,  we  can  get  about  150W  output  power,  but  it  will 
definitely  need  larger  volume  and  heavier  weight. 

3.3.  Effects  of  convection  heat  transfer  coefficient  of  each  side 

The  variation  of  performance  with  convection  heat  transfer 
coefficient  fq  or  h2  is  illustrated  in  Fig.  6.  Along  with  either  increase 
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Tc  (“o  Tc  CC) 


h i  or  h2,  the  output  power  and  efficiency  increases,  but  there  are 
some  differences.  When  h-,  or  h2  is  very  small,  both  the  change 
trends  of  output  power  and  efficiency  are  almost  the  same,  grow¬ 
ing  rapidly.  The  normal  triangular  curve  is  slightly  better  than  the 
circle  dot  curve  when  ft,  or  h2  is  lower  than  100  (W/m2I<). 
However,  as  the  convection  heat  transfer  coefficient  exceeds  100 
(W/m2  K),  it  begins  to  change  significantly.  The  growth  of  output 


power  and  efficiency  resultant  from  the  increment  of  hi  is  still  ra¬ 
pid,  while  the  growth  of  output  power  and  efficiency  resultant 
from  the  increment  of  h2  is  slow.  Although  the  temperature  differ¬ 
ence  grows,  the  temperature  of  the  hot  side  Ti  reduces  with  an 
increment  of  h2.  Considering  the  assumed  Seebeck  coefficient  de¬ 
creases  with  the  decrement  of  temperature,  so  the  improvements 
of  output  power  and  efficiency  may  be  partially  (even  all)  offset 
by  a  reduction  of  Ti.  On  the  contrary,  T,  increases  with  the  incre¬ 
ment  of  hi  and  the  temperature  difference  also  grows,  so  it  results 
in  significant  improvements  of  output  power  and  efficiency.  It  is 
noted  that  the  improvements  of  output  power  and  efficiency  will 
reduce  rapidly  if  the  temperature  of  the  hot  side  T,  is  higher  than 
the  best  working  temperature  range  of  the  material. 

In  fact,  it  is  much  easier  to  enhance  h2  than  hi,  such  as  using 
water-cooling  instead  of  air-cooling.  However,  not  enough 
improvement  can  be  achieved  by  only  increasing  h2. 

3.4.  Effect  of  height  of  the  PN  couple 

The  variation  of  performance  with  height  of  the  PN  couple  is 
shown  in  Fig.  7.  We  can  see  that  there  must  be  an  optimum  H  that 
maximizes  the  output  power.  When  other  parameters  do  not  vary, 
the  larger  the  H,  the  larger  the  temperature  difference  between  the 
two  sides  of  the  TEG;  however,  H  has  a  direct  ratio  relationship 
with  internal  resistance,  and  may  reduce  output  power  and 
efficiency. 


Fig.  6.  Effects  of  convection  heat  transfer  coefficient  of  either  s 
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Fig.  7.  Effect  of  height  of  PN  couple  on  output  power  and  Efficiency. 


The  dimensionless  figure  of  ZT  merit  of  semiconductor  material 
can  be  defined  as: 


where  T  is  the  Kelvin  temperature. 

There  are  three  ways  to  improve  ZT  value,  as  Eq.  (24)  shows,  and 
will  be  analyzed  respectively.  When  reducing  2,  the  best  output 
power  moves  toward  the  direction  of  the  reduction  of  H;  while 
increasing  the  S  or  a ,  the  best  output  power  moves  toward  to  the 
direction  of  the  increment  of  H,  and  those  optimum  output  powers 
are  nearly  equal.  All  efficiencies  increase  with  the  increment  of  H. 
Generally  speaking,  a  larger  Seebeck  coefficient  S  or  electric  con¬ 
ductivity  a  results  in  a  higher  optimum  H  value;  smaller  thermal 
conductivity  X  results  in  a  smaller  optimum  H  value.  Moreover, 
higher  height  of  material  usually  means  higher  cost  in  the  same 
cross-sectional  area  condition,  and  will  occupy  more  space  and 
weight.  Therefore,  when  we  design  a  TEG  device,  we  should  choose 
an  optimum  H  according  to  all  parameters  and  economy.  Just 
increasing  the  height  of  the  material  is  not  a  good  ideal  to  achieve 
greater  output  power  and  efficiency. 

3.5.  Effect  of  the  ratio  of  external  resistance  and  internal  resistance 

Fig.  8  shows  the  variation  of  performance  with  the  ratio  of 
external  resistance  to  internal  resistance.  A  maximum  output 
power  will  occur  at  the  ratio  of  about  1.5,  as  the  normal  triangular 
curve  shows.  This  is  different  from  a  common  electric  circuit  con¬ 
dition  when  external  resistance  is  equal  to  internal  resistance,  be¬ 
cause  the  generated  voltage  is  influenced  by  current.  Current 


decreases  when  resistance  increases,  causing  the  decrement  of 
heat  absorption  or  heat  release  of  the  Peltier  effect,  and  causing 
the  temperature  increment  difference  of  the  two  sides  of  the 
TEG,  leading  to  voltage  increases  according  to  the  Seebeck  effect. 
Same  analysis  is  done  for  the  ZT  value.  We  can  conclude  that  the 
inverted  triangle  improves  most  than  others,  and  the  most  likely 
reason  is  its  H  is  far  from  the  optimal  value,  so  income  is  greatest. 
When  H  is  near  optimum  value,  the  maximum  output  power  may 
be  almost  equal.  As  shown  in  Fig.  8,  we  can  conclude  when  we  dou¬ 
ble  ZT  value  in  the  three  ways,  best  output  power  and  efficiency  all 
move  toward  the  direction  of  the  ratio  increment,  and  the  reduc¬ 
tion  of  the  X  condition  moves  most.  Combined  with  3.4,  we  have 
observed  that  the  effect  of  S  or  a  is  almost  the  same,  while  having 
some  differences  with  the  effect  of  X  for  the  performance  of  the 
TEG. 


3.6.  Comparison  of  different  cooling-fluids 

Considering  that  the  convection  heat  transfer  coefficient  of  air  is 
smaller  than  that  of  water,  the  performance  of  the  TEG  using 
the  engine  cooling  water  as  cooling  fluid  is  investigated,  and  a 
comparison  analysis  is  carried  out  with  different  cooling  condition 
outcomes. 

The  vehicle  engine  cooling-water  temperature  is  generally 
between  80  °C  and  90  °C,  so  we  assumed  engine  cooling-water, 
temperature,  mass  flow  rate  and  coefficient  of  convective  heat 
transfer  respectively  as  80  °C,  200  g/s  and  1000W/(m2I<)  (coeffi¬ 
cient  of  the  convective  heat  transfer  of  water  is  greater  than  air’s 
in  uniform  conditions).  As  Fig.  9  shows,  the  temperature  of 
exhaust  gas  decreases  from  500  °C  to  340.6  °C,  otherwise  water 
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ratio  of  external  resistance  to  internal  resistance  on  output  power  and  efficiency. 


Fig.  8.  Effect  of  the 
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temperature  raises  from  80  °C  to  87.7  °C,  and  the  temperature  of 
semiconductor  material  ranges  from  213.4  to  166.5  °C  on  the  hot 
side,  and  from  107.3  to  104.9  °C  on  the  cold  side.  Output  power 
is  113.3  W  (improvement  of  72.7%  compared  to  air-cooling),  and 
efficiency  is  3.35%.  Water-cooling  is  obviously  better  than  air- 
cooling.  Using  the  vehicle  engine  water  coolant  as  the  TEG’s 
cooling-fluid  to  achieve  water-cooling  circulation  is  worthwhile 
in  practical  design  and  application. 


3.7.  Using  phase-change  material  for  heat  transfer  at  the  hot  side 

Combining  Figs.  3  and  9,  we  can  conclude  that  low  coefficients 
of  convective  heat  transfer  (especially  in  exhaust  pipe)  are  one  of 
the  important  limiting  factors  on  a  TEG  using  vehicles  exhaust  as 
resource.  Considering  that  backpressure  of  the  engine  must  affect 
its  efficiency,  it  is  quite  difficult  to  get  significant  improvement 
of  the  exhaust  coefficient  of  convective  heat  transfer.  As  phase- 
change  heat  transfer  has  a  huge  coefficient,  so  we  try  to  use 
phase-change  material  to  reform  our  experiment’s  design.  As 
Fig.  1 0  shows,  phase-change  fluid  boils  after  absorbing  heat  from 
the  gas  in  the  exhaust  pipe,  then  it  condenses  to  liquid  after  releas¬ 
ing  heat  to  the  hot  side  of  the  TEG.  The  condensate  liquid  flows  to 
the  boiling  side,  moving  cyclical.  Fig.  1 1  shows  the  comparison  of 
output  power  and  efficiency  between  a  system  using  phase-change 
material  and  a  normal  heat  convection  system  (we  assume  the 
phase-change  heat  transfer  coefficient  of  exhaust  is  5000 
W/(m2  K),  while  the  cooling  water  coefficient  of  convective  heat 
transfer  is  1000W/(m2K),  and  the  cooling-water’s  inlet  tempera¬ 
ture  is  80  °C,  and  mass  flow  rate  is  200  g/s). 

Using  material  with  lower  phase-change  temperature  (than  ex¬ 
haust  gas)  must  reduce  the  “quality”  of  exhaust  waste  heat,  but  a 
huge  heat  transfer  phase-change  coefficient  still  brings  significant 
improvement  of  both  output  power  and  efficiency,  as  Fig.  11 


Number  of  PN  couple  (j) 

Fig.  9.  Temperatures  of  each  unit  under  water-cooling  condition. 


boiling  exhaust  gas 


shows.  Circle  dot  is  the  calculation  result  using  phase-change 
material.  We  can  conclude  that  output  power  improves  greater 
than  10  times  at  45  g/s  mass  flow  rate  and  of  over  165%  at  5  g/s 
mass  flow  rate,  while  efficiency  improves  greater  than  148%  at 
45  g/s  mass  flow  rate,  and  of  over  313%  at  5  g/s  mass  flow  rate 
for  the  system  using  phase-change  material  versus  the  normal  sys¬ 
tem.  The  comparisons  of  exhaust  outlet  temperature  of  the  TEG 
with  different  mass  flow  rates  for  those  two  systems  are  shown 
in  Fig.  12.  We  observe  that  when  mass  flow  rates  are  below  35  g/s, 
the  exhaust  outlet  temperatures  of  the  TEG  are  equal  to  300  °C.  This 
means  the  heat  absorbed  by  the  phase-change  fluid  is  finite,  only 
part  of  PN  couples  are  working  (the  influence  between  PN  couples 
is  not  discussed  in  this  paper),  output  power  increases  greatly  with 
flow  rate  increases.  The  outlet  temperature  of  cooling-water  also  in¬ 
creases,  so  the  efficiency  decreases  slightly.  When  flow  rates  of  ex¬ 
haust  are  beyond  35  g/s,  the  output  power  of  the  TEG  reaches  a 
maximum  value,  and  increasing  the  flow  rate  does  not  have  influ¬ 
ence  on  output  power  and  efficiency,  while  the  exhaust  outlet  tem¬ 
perature  of  the  TEG  increases.  This  is  because,  limited  by  the 
dimension  of  the  TEG  device,  the  heat  absorption  capacity  of  the 
TEG  reaches  maximum,  all  PN  couples  are  working  within  large  tem¬ 
perature  differences. 

Without  question,  if  we  use  material  with  a  300  °C 
phase-change  temperature,  the  process  from  500  °C  to  300  °C  must 
reduce  the  “quality”  of  waste  exhaust  heat,  and  if  we  choose  mate¬ 
rial  with  a  higher  phase-change  temperature,  we  therefore  cannot 
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Fig.  10.  Simple  systematic  diagram  using  phase-change  i 
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Temperature  difference  (K) 

Fig.  13.  Comparison  of  present  simulation  result  and  experimental  result  ii 


recover  exhaust  water  heat  below  the  phase-change  temperature. 
It  will  not  be  a  small  loss.  So  we  put  forward  to  a  new  idea  using 
a  two-stage  (or  more)  combination  of  materials  with  different 
phase-change  temperatures,  along  the  flowing  direction  of  exhaust 
gas  for  cascade  utilization  of  exhaust  waste  heat  to  remit  this  prob¬ 
lem.  Detailed  analysis  of  this  ideal  will  be  reported  in  another 
paper. 

In  fact,  in  the  circumstance  using  phase-change  materials  to  im¬ 
prove  the  output  power  and  efficiency  of  the  TEG,  it  is  not  enough 
to  consider  just  the  temperature  of  the  phase-change  materials  and 
phase  change  heat  transfer  coefficient.  The  heat  convection  coeffi¬ 
cient  between  the  exhaust  gas  and  evaporator  is  also  an  important 
influencing  factor.  In  the  analysis  of  this  paper,  it  is  assumed  the 
heat  exchange  surface  between  the  exhaust  gas  and  evaporator  is 
large  enough,  compared  to  that  of  a  semiconductor  generator. 
So  the  influence  of  the  boiling  heat  transfer  coefficient  of 
phase-change  materials  and  the  convection  heat  transfer  coeffi¬ 
cient  between  the  exhaust  gas  and  the  evaporator  can  be  ignored. 


In  order  to  verify  the  validity  of  the  model  in  this  paper,  the 
present  simulation  result  and  the  experimental  result  in  literature 
[29]  are  compared,  as  shown  in  Fig.  13,  with  abscissa  for  tempera¬ 
ture  difference  and  ordinate  as  output  power.  As  can  be  seen  in  the 
temperature  range  shown  in  the  figure,  good  agreement  exists 
with  those  obtained  in  this  paper,  so  the  model  calculation  results 
are  reliable. 

4.  Conclusion 

Based  on  Fourier’s  law  and  the  Seebeck  effect,  this  paper  pre¬ 
sents  a  mathematical  model  of  a  Thermoelectric  Generator  (TEG) 
device  using  the  exhaust  gas  of  vehicles  as  heat  source,  and  preli¬ 
minary  analysis  of  the  impact  of  relevant  factors  on  the  output 
power  and  efficiency  of  the  TEG.  Main  conclusions  are  as  follows: 

(1)  In  an  air-cooling  condition,  the  increasing  convection  heat 
transfer  coefficient  of  either  side  of  the  TEG  will  result  in 
output  power  and  efficiency  increases,  but  the  effects  are 
different.  The  change  trends  of  the  output  power  and 
efficiency  are  roughly  the  same  when  the  heat  transfer  coef¬ 
ficient  is  small.  When  either  convection  heat  transfer  coeffi¬ 
cient  exceeds  the  other,  the  improvements  of  output  power 
and  efficiency  by  the  increment  of  convection  heat  transfer 


coefficient  of  the  high-temperature-side  are  greater  than 
by  the  increment  of  convection  heat  transfer  coefficient  of 
the  low-temperature-side.  But,  because  of  the  limited  space 
and  backpressure  of  the  exhaust,  the  increment  of  convec¬ 
tion  heat  transfer  coefficient  of  the  high-temperature-side 
is  quite  difficult  or  costly. 

(2)  With  the  change  of  height  of  the  PN  couple,  there  is  a  peak 
value  of  output  power.  This  best  height  value  for  peak  point 
decreases  with  the  decrement  of  thermal  conductivity  of 
material,  and  increases  with  the  increment  of  the  Seebeck 
coefficient  and  electric  conductivity  of  the  material. 

(3)  The  output  power  and  efficiency  of  the  TEG  change  with  the 
change  of  the  ratio  of  external  resistance  to  internal  resis¬ 
tance,  and  there  also  exists  a  maximum.  The  maximum 
value  appears  when  the  external  resistance  is  greater  than 
the  internal  resistance.  It  is  different  from  a  common  circuit 
at  the  time  the  external  resistance  is  equal  to  internal  resis¬ 
tance.  And  with  the  increment  of  the  ZT  value,  the  maximum 
value  ratio  moves  toward  the  direction  of  increasing  the 
ratio  of  external  resistance  to  internal  resistance. 

(4)  A  preliminary  analysis  of  the  possibility  of  using  phase- 
change  materials  to  reinforce  the  heat  transfer  property  of 
the  high-temperature-side  of  the  TEG  was  carried  out,  and 
the  results  show  that  using  appropriate  phase-change 
materials  can  significantly  improve  the  output  power  and 
efficiency  of  the  TEG. 
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